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Interpretation of body-mounted 
accelerometry in flying animals and 
estimation of biomechanical power 

R. J. Spivey and C. M. Bishop 

School of Biological Sciences, Bangor University, Bangor, Gwynedd LL57 2UW, UK 

An idealized energy fluctuation model of a bird's body undergoing horizon- 
tal flapping flight is developed, focusing on the biomechanical power 
discernible to a body-mounted accelerometer. Expressions for flight body 
power constructed from root mean square dynamic body accelerations and 
wingstroke frequency are derived from first principles and presented in 
dimensionally appropriate units. As wingstroke frequency increases, the 
model generally predicts a gradual transition in power from a linear to an 
asymptotically cubic relationship. However, the onset of this transition 
and the degree to which this occurs depends upon whether and how for- 
ward vibrations are exploited for temporary energy storage and retrieval. 
While this may vary considerably between species and individual birds, it 
is found that a quadrature phase arrangement is generally advantageous 
during level flight. Gravity-aligned vertical acceleration always enters into 
the calculation of body power, but, whenever forward acceleration becomes 
relevant, its contribution is subtractive. Several novel kinematic measures 
descriptive of flapping flight are postulated, offering fresh insights into the 
processes involved in airborne locomotion. The limitations of the model 
are briefly discussed, and departures from its predictions during ascending 
and descending flight evaluated. These findings highlight how body- 
mounted accelerometers can offer a valuable, insightful and non-invasive 
technique for investigating the flight of free-ranging birds and bats. 



1. Introduction 

Birds flap their wings in order to achieve weight support and locomotion [1-3]. 
Experiments using high-frame-rate video footage to monitor wing and body 
motions of birds or bats flying in wind tunnels have combined the findings with 
aerodynamic results and accelerometry to estimate overall energy expenditure 
during flight [4-6]. Such approaches have been experimentally valuable and theor- 
etically illuminating, enabling the refinement of aerodynamic theory, but in studies 
involving free-ranging animals where trailing wires and heavy equipment cannot 
be used, ambulatory recording of body acceleration offers a viable and practical 
alternative. Accelerometry was initially restricted to wind tunnel work [7], 
but has now been miniaturized and demands relatively little electrical power. 
Commercially available micro-electromechanical transducers are now capable of 
faithfully recording high-frequency vibrations, offering a new means of studying 
the characteristics, kinematics and energetics of free-ranging avian flight [8,9] 
and, indeed, animal locomotion in general [10]. Because the long-term study of 
birds in the wild is becoming increasingly feasible, there is new scope to assess 
some of the difficult choices birds face during long-range migrations [11,12]. 

Traditional techniques for monitoring the metabolic rate or power input of free- 
ranging vertebrates include doubly labelled water [13,14] and heart rate,/ h , derived 
from electrocardiography [15,16]. The latter approach offers good temporal resol- 
ution but has historically necessitated calibration of / h against measurements of 
oxygen consumption, Vq 2 • However, the direct translation of / h to Vq 2 may now 
be possible for endotherms undergoing primary mode locomotion if augmented 
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by knowledge of heart and body mass [17]. Accelerometry has a 
similar potential to monitor instantaneous biomechanical power 
output [18] during locomotion, and strong correlative relation- 
ships between body acceleration and Vq 2 have been found in 
animals running on treadmills [18,19]. As one might expect, 
body accelerations during flight are generally elevated compared 
with other forms of locomotion [14,20-22]; however, a theoreti- 
cal understanding of how body acceleration relates to the 
biomechanical power of flapping flight has not yet been eluci- 
dated. This study aims to address this by setting out a 
mathematical model that assists the interpretation of accelero- 
metry data captured from birds during flight. Novel measures 
descriptive of flight kinematics, integral to this modelling, shall 
also be derived. 

Instruments that log acceleration are, for practical reasons, 
generally attached to the torso of a flying animal. As wings 
are coupled to the body, this offers much promise as a non- 
invasive tool that can help estimate the biomechanical 
power (and indirectly or proportionally, the metabolic 
costs) associated with flight [18,19,23], with the potential to 
augment or replace existing methods [15,24]. Activity-related 
accelerations can be decomposed into the sum of dynamic 
and static accelerations which can be respectively derived 
by high- and low-pass filtering either in the time or the fre- 
quency domain. To date, biologists have found the dynamic 
component most informative with regard to correlations 
with energy expenditure. Two time-averaged measures of 
dynamic body acceleration (DBA) have been used when 
studying the energetics of animals [18,25]. Overall dynamic 
body acceleration (ODBA) is a running average of the 
L 1 -norm of the dynamic acceleration [10,18]. The variant 
using the L 2 -norm, which accurately encapsulates vectorial 
length, is known as vectorial dynamic body acceleration 
(VeDBA) [26,27]. An immediate difficulty with using any 
acceleration measure as a proxy for estimating biomechanical 
power in the absence of empirical calibration is that the funda- 
mental units of acceleration, namely IT -2 , are different both to 
those of power, ML 2 T~ 3 , and those of mass-specific power, 
L 2 T~ 3 . Overcoming this inevitably requires the development 
of some theoretical model descriptive of the biomechanics of 
bird flight with respect to body acceleration. 

Under the hypothesis that DBA is closely related to over- 
all metabolic costs, accelerometry has been successfully 
applied to a wide variety of other animals [10]. Consequently, 
it is not unreasonable to expect that a correlation could also 
exist between DBA and the biomechanical power directly dis- 
cernible using a body-mounted accelerometer (body power) 
during steady horizontal flight and the metabolic rates of 
birds during flight (subject to the additional uncertainties 
of the value for the mechanochemical conversion efficiency 
of muscle [28]). Therefore, the ansatz is adopted here that 
the kinetics of the body should reflect the kinetics of the 
wings, thereby ultimately allowing the biomechanical costs 
incurred during avian locomotion to be estimated. This 
work primarily focuses on the interpretation of data from 
accelerometers attached to the body of a flying bird in the 
absence of additional information, a constraint demanding 
the development of some mathematical model to theoreti- 
cally bridge the divide between body vibrations and overall 
biomechanical costs. The relationship between decomposed 
vertical and horizontal dynamic accelerations, and the var- 
ious components of the energy associated with the body are 
investigated, lateral components being neglected due to 



the symmetrical beating of the wings. Birds must find ways 
of contending with the weight of gravity when airborne 
and flapping flight demands significant energy expenditure 
[1,3]. Motions of the body on the vertical axis differ from 
motions within the horizontal plane as they involve changes 
in gravitational potential. Thus, if the biomechanical power 
during flight is to be estimated solely from accelerometry, 
then it is essential to pay heed to the direction of gravity 
and hence also the absolute orientation of the transducer. 
The ultimate aim of this work is to derive estimates of 
acceleration-based biomechanical body power during flap- 
ping flight from first principles, potentially leading to the 
future interpretation of accelerometry from flying animals 
without the need for direct calibration. Novel statistical 
measures derived here may also be informative of flight 
kinematics, pertaining to energetically significant transitions 
in the wingstroke frequency, such as the relative phase and 
amplitude of forward and vertical body oscillations. 

2. Developing a model 
2.1. Preliminaries 

Due to the pulsatile character of avian flight associated with 
the periodic contraction of powerful wing muscles, the energy 
associated with each wingstroke is delivered sporadically. 
Efficient flight confers evolutionary advantages, so the effort 
required by a bird to flap its wings is likely to achieve useful 
goals such as forward propulsion against aerodynamic drag, 
the countering of gravity, changes in velocity, ascent/ descent 
and general manoeuvring. The mechanical energy of a bird 
will fluctuate in time and in still air there is a metabolic cost 
when the total energy increases. When it decreases, the dissipa- 
tion of mechanical energy into the surrounding air is used to 
accomplish these various flight goals [29]. In the absence of grav- 
ity and an atmospheric medium, the mean mechanical power 
required by a vibrating but dissipationless mechanical system 
is precisely zero, because the total energy of the system is 
constant at all times. However, energy will necessarily be trans- 
formed or exchanged between different elements of the system 
within individual vibration cycles. Birds, on the other hand, 
remain airborne and sustain forward momentum despite air 
resistance by doing mechanical work which they never recover. 
Notwithstanding this, birds may be able to temporarily store 
and retrieve energy within individual wingbeat cycles in a simi- 
lar manner to an idealized dissipationless system, so this should 
be considered when developing the model, as such a mechanism 
may provide significant flexibility to execute flight more 
efficiently under certain circumstances. 

We assume a triaxial accelerometer of negligible mass 
securely affixed to the body of a bird undergoing horizontal 
flight at a steady air speed. It has been known for several 
decades that, provided adequate consideration is given to har- 
ness design, accelerometers can be externally mounted to birds 
with negligible oscillation relative to the body [30]. Raw data 
from an accelerometer can be reoriented using mathematical 
transformations, an issue returned to later, so this analysis 
proceeds by taking the z-axis to be oppositely aligned to gravity 
and the y-axis to correspond to the direction of forward motion. 
An inertial frame of constant velocity comoving with the bird is 
considered. The majority of the power in the measured accelera- 
tions resides at the fundamental wingstroke frequency /, the 
angular frequency of the wingstrokes being o) = 27rf. This is a 
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Figure 1. A typical power spectrum derived from the Fourier decomposition 
of about 10 s of vertical axis accelerometry data captured by the authors at 
400 s s -1 from the body of a freely flying homing pigeon during horizontal, 
straightline flight while returning to a loft. The peak at the fundamental 
wingstroke frequency (at around 6.5 Hz) dominates over that of higher 
harmonics, suggesting that the acceleration on the bird's body is not only 
periodic but essentially sinusoidal. 

crucial assumption, backed up by the acceleration power 
spectrum obtained via accelerometry experimentally collected 
from a flying pigeon, as presented in figure 1. Because power 
is concentrated at the fundamental frequency, the oscillations 
of the body in the vertical axis can be accurately approximated 
by a sinusoid with maximum excursion B relative to the mean 
altitude. The vertical displacement z and vertical velocity z 
therefore vary as follows: 



and 



B cos(cot) 



-Bcosm(o)t). 



(2.1) 
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where c = yja 2 + b 2 — lab cos 4> and H is the Heaviside step 
function in which H( £) = 1 if £ > 0, otherwise H( £) = 0. It 
can be seen that the resultant sinusoid will generally have 
an intermediate phase shift. 

2.2. Relative phase between vertical and horizontal 
oscillations of the body 

Vibrations on both axes will exhibit simple harmonic motion 
with both K y and K z varying at the same frequency. However, 
before an expression for the variation in y can be written as a 
sinusoid, there is a need to carefully consider its relative 
phase 0 with respect to z. Even if the bird has considerable 
freedom to adjust the phase of its motion on the horizontal 
axis, arrangements that minimize the power required in 
order to sustain flight confer evolutionary advantages. Intro- 
ducing a horizontal vibration amplitude A, the counterpart to 
B for the vertical axis, one can write 



and 



y = Acos(cot - 



y = —A(osm(cot + 6). 



(2.4) 
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The convention is adopted that neither A nor B can be nega- 
tive. The total energy associated with the bird's body, E h (t), is 

- b y 2 



m^z m^y 



E h (t) = U + K z + K y = m h gz + — — + 
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and 
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It is assumed that no energy is stored either elastically or 
as rotational kinetic energy, leaving the gravitational poten- 
tial energy U = m h gz, the kinetic energy associated with the 
vertical axis K z = v%z 2 /2, and the kinetic energy associated 
with the forward axis K y = m^y 2 /2 to be considered. Here, 
m h is the mass of the body of the bird, and g is the gravita- 
tional acceleration. The beating of both wings is assumed to 
always be in phase with one another. Because of the bilateral 
symmetry in birds, we neglect any kinetic energy associated 
with the x-axis. 

The wings also have gravitational and kinetic energies 
(both vertical and horizontal). However, these are not 
directly estimable from a body-mounted accelerometer, so 
terms for these quantities are not explicitly included in the 
analysis. Instead, it is assumed that the activity detected by 
the accelerometer will, to a first approximation, be an attenu- 
ated reflection of the total biomechanical output of the bird. It 
is reasonable to think that energy associated with the wings 
can be subsumed into existing terms by constructive or 
destructive superposition. This follows from the fact that 
the addition of two arbitrarily scaled sinusoids of the same 
frequency, but different phase results in a rescaled sinusoid 
of the same frequency: 



a sin(orf) + b sin(cot + cp) = c sin 



tot + arctan 



/ b sin 4> 
\a + b cos </> 



+ 7rH(-a — bcos <j>) 



(2.3) 



During each wingstroke, the flight muscles of the bird 
must supply energy when E h is increasing. An important 
aspect of the present model is the assumption that when E h 
decreases, energy is irrecoverably lost to the environment 
and exploited so as to achieve propulsion and weight 
support. This can be compared with aerodynamic models 
which assume the kinetic energy of the wings is never 
recovered [31,32]. The mean power associated with the 
body, (P b ), can then be determined according to 

(Pb) = E h (t)H(E b (t)) dt (2.8) 



K y = 0 and 



Consider, for now, the case in which A = 0. One then has 
E h = U + K z . Gravitational potential energy 
attains a minimum at t = it I co when z = — B and a maximum 
at t = 0 when z = B. At both these times, K z reaches its mini- 
mum of zero. Because sin 2 ^ = (1 — coslcot)/!, the vertical 
kinetic energy varies sinusoidally at twice the frequency of 
the gravitational energy. Thus, E h (t = it I co) is always a mini- 
mum. Although U is always a maximum at t = 0, it is 
possible that E h is not a maximum at that time if K z thereafter 
increases more rapidly than U decreases. This possibility is 
apparent upon inspection of the time derivatives of E h/ the 
roots of which correspond to stationary points: 

Eb(OL=o = m\jB(osm(cot)[B(o 2 cos(cot) — g] (2.9) 



and 



EbWL=o = m h B(o 2 [Bco 2 cos(2cot) -gcos(cot)}. (2.10) 



Owing to the sm(cot) term, stationary points exist at t = 0 
and t = it I oo. At t = it/ oo, one has Eb = nh 0 Boo 2 {Boo 2 + g 2 ) > 0, 
so this stationary point is always a minimum. However, 
E\y = m\yBcx) 2 (Bco 2 — g) for the stationary point at t = 0, 
which corresponds to a maximum when co 2 < g/B and a 
minimum when co 2 >g/B. If both these are minima, 
additional stationary points occur when cos(cot) = g/Boo 2 . 
Using cos(2cot) = lg 2 /B 2 oo 4 — 1, one then finds that at those 
times Eb = ^b(g 2 — B 2 oo A ) < 0, confirming that these points 
correspond to maxima of E b . 

Therefore, when co 2 < g/B, the mean body power is simply 
<P b ) = ((o/Itt) (Eg** - E^), where E™* = E b (0) = m h gB and 
Eg™ 1 = E^(tt/oo) = — nbogB. However, when oo> g/B, then 
due to the existence of a new maximum in E h at t ^ 0, E b nax 
will exceed E b (0), and the mean power will inevitably rise. 
During intense flight, there will be a high wingstroke frequency, 
and for each wingstroke the variation in kinetic energy will 
increase and the variation in gravitational energy will decrease, 
allowing the variations in kinetic energy to become dominant. 
However, if it were possible to temporarily store some of the ver- 
tical kinetic energy as horizontal kinetic energy, and retrieve it 
later in the wingstroke cycle, then this elevated maximum in 
E h could be avoided, and the concomitant increase in power 
is eliminated. 

The transition between the two regimes occurs at g = Boo 2 , 
corresponding to the peak gravitational energy U max = m h gB 
being equal to twice the peak vertical kinetic energy 
2K™ X — m\)B 2 (x) 2 . Clearly, if the total kinetic energy K yz = 
K y + K z did not fluctuate at all, then K z could be arbitrarily 
large without incurring any additional increase in mean 
power. The total kinetic energy is 
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'A 2 + B 2 -A 2 cos(loot + 20) - B 2 cos(loot)]. 
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If K yz is constant, then its derivative 



Ky Z — K Z + Ky 



^ [B 2 sin 2 (cot) + A 2 sin 2 1 



cot 



and 



nib co 



[B 2 + (A 2 -B 2 )cos 2 (wf)] 



tJB 2 



(2.16) 



(2.17) 



A relative phase shift of ±tt/1 between the horizontal and 
vertical axes corresponds to what is commonly termed a quad- 
rature phase arrangement. Only in this circumstance does it 
hold that Ky is a maximum when K z is a minimum, and vice 
versa. This maximizes the potential for shuttling energy back 
and forth between the two axes, a useful property that the 
bird might be able to exploit to decrease its mean power. 
Although K y will be maximal at t= 0 and t = it/ oo, because 
these maxima are equal, they have no effect on the difference 
E^(t — 0) — E h (t = it/ co). In fact, this holds for any value of 0 
because sin 2 x = sin 2 (x + it), but K yz is constant only when K y 
and K z are in antiphase. As the kinetic energies vary at 
double the fundamental wingstroke frequency, this occurs 
when 6 = ±it/2. 

While variations in K yz can be completely eliminated, it 
may not be necessary for the bird to do so because, as will 
be shown, the bird can in some cases also minimize power 
when A/B < 1 with excessive fluctuations in K z being com- 
pletely tamed by smaller fluctuations in K y . This may be 
preferable as it reduces extraneous energy losses and amelio- 
rates the vibrations transmitted to the head of the bird, which 
might otherwise make flight an unnecessarily uncomfortable 
experience compromising visual acuity [33]. Note also that if 
A/B > 1, then fluctuations in K y may not be adequately 
absorbed by fluctuations in K z . 

Because birds are expected to have a maximum forward 
velocity at the end of the downbeat, the phase that makes 
y(0) maximal is chosen (0= —tt/2). Hence, 6 can be 
eliminated from the expressions for y and y by writing 

y = AsinO£) (2.18) 

and 
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[A 2 sin(2cot + 20) + B 2 sin(2orf) 



(2.13) 



must be zero at all times. It is apparent from this expression 
and (2.3) that K yz varies sinusoidally. The amplitude of these 
fluctuations vanishes when 

A 2 sin loot cos 20 + A 2 cos loot sin 20 + B 2 sin loot = 0 (2.14) 

and 



y = Aoocos(cot). 



(2.19) 



The upbeat commences at t = 0 when y is maximal. The 
body is then at its maximum height above the ground. 
The body and the wings are in antiphase on the vertical 
axis, and the same should also be approximately true of the 
forward axis. To summarize the findings of this section, 
power reduction is possible only when g/Boo 2 < 1 and is 
best achieved by a quadrature phase arrangement. 



tan 2oot - 



sin loot 



-A 2 sin 20 



cos2oot A 2 cos 20 + B 2 ' 



(2.15) 



Because 0 is constant, this would demand that t is also 
constant. However, K yz can be zero, if the sin2co£ and the 
cos2oot terms are simultaneously zero. One then finds that 
the cos2oot term vanishes if sin 20 = 0, which is satisfied 
when 0 = nir/2 (n G Z). The sm2oot term is zero when 
A 2 cos 20 + B 2 = 0, which yields 0 = arccos(-B 2 /A 2 )/2. When 
0= nir, one has cos 20=1 which must be rejected as it pre- 
dicts A 2 =-B 2 . However, 0=(n± 1/2)tt yields A 2 = B 2 , 
which is acceptable. Therefore, for K yz to remain constant 
requires A = B and 0=±tt/1. This can be seen in the 



2.3. Quadrature phase flight 

Unless expressly stated otherwise, the analysis now proceeds 
by assuming quadrature phase applies. In order to assess the 
biomechanical power in the body, one is interested in deter- 
mining the maxima and minima (stationary points) of the 
body energy, necessitating looking for roots of the first-time 
derivative and inspecting their signs by taking the second- 
time derivative. Both K yz and U vary sinusoidally but because 
they do not vary at the same frequency their sum is not a 
simple sinusoid, demanding that calculus be used. Body 



energy now reads 

EbW = (^)[2gBcos(ttrf) + B 2 o) 2 sin 2 (cot) + A 2 o? cos 2 (cot)] , 

(2.20) 

and the first derivative is 

E b = m\>o) 3 (B 2 — A 2 ) sin(ftrf) cos(cof) — m b gB(w sin (cof) (2.21) 
and 

E b = m h cosm(cot)[oo 2 (B 2 - A 2 ) cos(cot) - gB]. (2.22) 

Stationary points exist when either cos(cot) = gB/ co 2 (B 2 — A 2 ) 
or sm(cot) = 0. The first condition has real solutions only if 
(o 2 > gB/\B 2 - A 2 \. When oo 2 > g/B, the smallest value of A 
that prevents a maximum in E b from arising anywhere but 
at t = 0 occurs when cos(cot) = 1, from which the smallest 
value of A that minimizes the power can be determined. If 
this is an overriding consideration with the need to maintain 
visual acuity a secondary concern, this value of A would 
seem to be optimal: 



optimal acceleration partitioning 
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opt 



BJl-g/Ba 



(2.23) 



Note that A opt must be real (because co 2 > g/B) and 
because A 2 pt < B 2 at all times, birds have no need to fly 
with A > B. Moreover, for given values of A and B, the opti- 
mum angular frequency comes to co opt = yJgB/(B 2 — A 2 ) 
which exceeds the critical frequency co c = yJg/B above which 
power can be reduced if A ^ 0. The optimum ratio of A/B 
ensures fluctuations in K y are sufficiently large to avoid 
excessive fluctuations in K z , its value depending only on 
the ratio co/ co c : 




Figure 2. The optimum value of A/B, the ratio of the horizontal and vertical 
vibration amplitudes, at which body power is minimized and dynamic acceleration 
on the forward axis attains its smallest acceptable value, see equation (2.24). 



Note that the stationary point in E b at t = it/ co has second 
derivative 



E h (7T/co) = m h oo 2 [(B 2 - A 2 )co 2 + 



(2.28) 



This is negative if gB < co 2 (A 2 - B 2 ) so that when A 2 > B 2 , 
which gives rise to stationary points at cos(cot) = 



co 2 (A 2 — B 2 ), the stationary point at t = it/ co becomes a maxi- 
mum. This is a point of inflection when gB = co 2 (A 2 — B 2 ) 
corresponding to the maximum value of A at which power is 
minimized: 
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This is plotted in figure 2. The second derivative of the 
body energy is 

E b ra b co 4 (B 2 - A 2 ) cos(2oot) - m h gBoo 2 cos(oot). (2.25) 



The condition that Eb — 0 at t = 0 yields A = A { 



opt/ 



showing that the maximum at t = 0 is then also a stationary 
point of inflection. Now, consider the stationary points that 
arise when cos(cot) = gB/ oo 2 (B 2 — A 2 ). Using the identity 
cos(2cot) = 2cos 2 (cot) — 1, one finds that 



E h =2m h co 4 (B 2 - A 2 ) 



? 2 B 2 



oo\B 2 -A 2 ) z 
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If A 2 < B 2 , the sign of Eb must be negative if gB < co 2 
(B 2 — A 2 ). This satisfies the condition required for the stationary 
points to exist, establishing that those corresponding to cos(cot) = 
gB/co 2 {B 2 - A 2 ) < 1 must be maxima if A 2 < B 2 . Considering 
the alternative situation in which A 2 > B 2 , the sign of Eb must 
be positive if gB < co 2 (A 2 — B 2 ). Hence, the stationary points at 
which cos(cot) = gB/ co 2 (B 2 - A 2 ) < 1 must then correspond to 



This has a similar form to the earlier expression 
A opt = co 2 — col/ oo. ki r d has minimized its mean body 
power if A opt < A < A max or equivalently if co 2 < gB/ |B 2 - A 2 \. 
When this condition is satisfied and co 2 > g/ B, flapping flight 
demands less power than when forward axis vibrations are 
absent altogether (A = 0). For a fixed value of B, this is illustrated 
in figure 3. 



3. Estimation of body power 

The calculation of mean body power here involves the inte- 
gration of E b only as it increases according to (2.8). This can 
be easily accomplished by subtracting the minima of E b 
from successive maxima of E b . If t = 0 at t = t 0 , t = tt/ co at 
t = h and t= ±axccos[gB/co 2 {B 2 - A 2 )\/co at t = U then 
only these times need be considered in order to determine 
the minima and maxima of the body energy variation. With 
reference to figure 4, it can be seen that there are three 
distinct cases to consider with stationary points located at 
the following times: 



case l:t = t 0 and t = h: {gB > oo 2 \B 2 -A 2 \}, 
case 2: t = t 0 , t = h and t = f* {gB < co 2 {B 2 - 

E h (Q > E h (t 0 )} and 
case 3: t = t 0 , t = h and t = U {gB < oo 2 (A 2 - 

E h (Q < E h (h)}. 
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Figure 3. Body power is minimized providing the ratio AIB falls within 
the interval (By/l - oJ-JoZ.B^jA + to]/ to 2 ) if co > co c or 
(0, B^/l + co\/co 2 ) if co < co c . The upper and lower limits to AIB 
each approach unity as colco c — > oo (see text). 
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Figure 4. Variation of body energy assuming motions on the forward and ver- 
tical axes are in quadrature phase (see text). In this model, the bird provides the 
energy required for the traces to increase with time but does not recover that 
energy when it declines. Thus, differences in energy between the stationary 
points (marked by dots) allow the calculation of mean power. There are two 
limiting cases: (i) co 2 = gBI(B 2 - A 2 ), when B 2 > A 2 , and (ii) co 2 = gBI 
{A 2 - B 2 ), when A 2 > B 2 (both plotted) delineating three modes of flight. 
The case A = B bisects the range in which the mean power attains a minimum, 
which becomes narrower as wingstroke frequency increases (figure 3). Outside 
this range, stationary points appear at times intermediate between the mid- 
point and the start/end points of a wingstroke cycle (at t = f*, upper and 
lower traces provide examples). (Online version in colour.) 



The body energy at t 0 and t\ is simply 



and 



Eb(fo)=(f)04V + 2gB) 



(3.1) 



(3.2) 



When t = U, sin (cot*) is obtained as 1 - cos (cot*): 
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After some algebra, the total energy of the body when 
f = f* is found to be 
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3.1. Case 1 

This case corresponds to the shaded region in figure 3 and the 
central traces of figure 4 for which the calculation of mean 
power is particularly simple: 
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3.2. Case 2 

This case corresponds to the zone below the shaded region 
in figure 3 (see also the uppermost trace in figure 4). The 
stationary point at t 0 is now a minimum and two new 
maxima arise at t = t*. Recalling that in this case, B 2 >A 2 . 
the power is 
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3.3. Case 3 

This case corresponds to the zone above the shaded region in 
figure 3 (see also the lowermost trace in figure 4). The station- 
ary point at t\ is now a maximum and new minima arise at 
t = U. One has A 2 > B 2 , and the power now comes to 

(Pbs) = ^ [E b (fo) - E„(f*) + E h (h) - £„(*♦)] P-S) 
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The poles appearing in the expressions for (Phi) and (Pb3) 
as A — > B are avoided because when co 2 > gB/\B 2 — A 2 \ one 
has A 2 ^ B 2 . It is apparent that the results for cases 2 and 3 
are equivalent but for the reversal of sign in the B 2 - A 2 
terms. Therefore, it would be acceptable to take the modulus 
of either expression without expressly checking whether 
A 2 > B 2 . In the limiting cases where gB = co 2 \B 2 - A 2 \, the 
prediction of case 1 coincides with that of case 2 or 3. For 
instance, when A = A opt = B^l -g/Bco 2 then B 2 - A 2 = 
gB/co 2 and it can be seen that (Phi) reduces to (Pbi)- 
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The importance of correctly distinguishing between 
case 1 and cases 2 and 3 is stressed, because (Phi) an d 
(Pb3) overpredict the true power when gB > co 2 \B 2 — A 2 \ 
and (P bl ) underpredicts the true power when gB < 



co z \B z 
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3.4. Deriving body power using data from 
accelerometers 

The root mean square (or r.m.s.) value of a quantity is a 
measure commonly used in physics and engineering and 
can lend itself to the description of DBA (table 1). As such, 
it offers an alternative to ODBA and VeDBA. The r.m.s. 
value of a discr ete set o f N dynamic single axis accelerations 



Table 1. Variables used in this study. 



a dyn is a rms — y #~[ yn /N, and that of a simple sinusoid such as 
\p = p sin(f ) is 



*Arms = V <^> 
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(3.11) 



Static and dynamic accelerations can be respectively 
derived from raw acceleration data using low-pass and 
high-pass filtering techniques. While the static acceleration 
is useful in determining the vertical, gravity-aligned axis, 
there are many circumstances where the direction of the 
forward axis is more ambiguous. However, birds will gen- 
erally adjust the roll of their bodies during flight so that the 
static acceleration vector remains dorsally aligned. It is 
therefore likely that the r.m.s. value of the dynamic sway, 
^rms/ will be appreciably smaller than the r.m.s. value of 
the dynamic surge, i/ rms . In principle, this allows for reor- 
ientation of accelerometry data during post-processing by 
application of a rotation matrix whose components can be 
inferred by analysis of the data. A method for achieving 
reorientation is outlined in appendix A. Estimates of the 
dynamic surge, i/ rms , and the dynamic heave, z rms , are 
readily obtained after high-pass filtering of the reoriented 
acceleration components. These relate to to, A and B as 
follows: 
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Hence, body power can be directly evaluated from r.m.s. 
heave and surge: 
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(3.15) 



These expressions conveniently obviate the need to 
double integrate acceleration data in order to obtain the 
values of A and B directly, which is generally challenging 
due to the baseline drift introduced when integrating. 
However, A and B each feature in the true-or-false test 
gB < oj 2 \B 2 — A 2 \ that determines which expression for 
power is valid. Therefore, a reformulation of the discrimi- 
nant is also desirable, and because i/ rms /z r ms = A/B it 
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follows that 
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It is therefore possible to determine, without knowledge 
of either A or B, nor indeed co, the regime in which the 
bird is flying. (Pbi) should be used when g'z rms > y/2 
l4ns -vLb\> otherwise (P b2/3 ) is applicable. 

The critical frequency ratio, co/ co c , can similarly be trans- 
lated using (3.13) and the knowledge that co c = y/g/B 



(3.17) 



It is also possible to express co/co opt using only r.m.s. 
accelerations: 
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Figure 5. Above a threshold frequency that depends on the ratio AIB, when 
A ^ B an asymptotically cubic response is inevitable at high frequencies. 
However, body power is always minimized when A = B. See equations 
(3.5), (3.7), (3.9) and (3.22). (online version in colour.) 



3.5. Relative body power 

According to this model, body power retains linearity with 
wingstroke frequency until, and as illustrated in figure 5, 
the angular frequency exceeds 
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This can also be conveniently gauged from r.m.s. accelera- 
tions alone because 



V2I 



■2/Ll 



(3.20) 



The expressions for (P b2 ) and (P^) can be recast as 
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(3.22) 



It is now apparent that when co > co 0/ the response 
becomes asymptotically cubic at higher wingstroke fre- 
quencies, i.e. (Pb2,3) oc <^ 3 . However, at lower wingstroke 
frequencies, flight proceeds within the linear regime of 
(3.5). For case 1, it is possible to write (Pbi) — mg 2 co/ ttco 2 , 
allowing the ratio of (Pbi)/ (^2,3) to be expressed in a particularly 
simple form 
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(3.23) 



The response is plotted in figure 6 and can be compared 
with figure 5 where the wingstroke frequency is expressed 
in units of co/ co c . If co < co 0 , then the expression for (Pb2,3) is 
not physically meaningful, and the mean power is always 
given by (Pbi). When co > co 0 , power can always be reduced, 
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Figure 6. The ratio (P^/iP^) as a function of normalized wingstroke 
frequency oj/ojq. When co < oj 0i forward vibrations are incapable of 
reducing the mean body power. Note that (P bl ) never exceeds (^2,3)- 

because the ratio (Pb2,3)/(Pbi) then exceeds unity. Further- 
more, it grows without limit as co — > 00. 

Power is always minimized when A = B but it is interesting 
to know how the ratio (Pbi^)/ (Pbi) grows for other values of A/ 
B when co > co 0 . First, note that co 0 is related to co c according to 
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(3.24) 



The relative flight cost ratio (Pbi^)/ (Pbi) can be expressed 
either in terms of co/ co 0 or A/B and co/ co c : 
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Figure 7. The relative cost of flapping flight normalized to that possible 
when A = B. For example, at a constant wingstroke frequency a)/a) c = 2, 
flight costs more than double when A/B < \^Jl — ^3 or 
A/B > \\/6 + y/l. (Online version in colour.) 



and 

aw = i 

(Phi) 2 
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A colour-coded contour plot of this function is presented 
in figure 7 (online). 



4. Non-quadrature phase 



We reiterate that this modelling pertains to steady, horizontal 
flight. The possibility exists that departures from quadrature 
phase may be advantageous during ascent or descent, but 
the mathematics in such cases is considerably more involved. 
Nevertheless, these situations were numerically investigated 
under the assumption that A = B with an additional power 
component representing the change in gravitational energy 
with time, whose mean value can be estimated using an alti- 
meter or GPS device [34,35]. The results are presented in 
figure 8. It can be seen that predictions of the quadrature 
phase model are still accurately upheld in most circumstances 
involving realistic rates of ascent or descent. 

Although the optimal arrangement is never far from 
quadrature phase, the error is sensitive to the rate of ascent, 
and is asymmetrical in that it grows faster with descent 
than ascent. This is evident from the plot presented in 
figure 9. It has been reported that the power requirements 
of moderate ascending and descending flight in pigeons 
can be accurately estimated by summing the power required 
for level flight with the rate of change of gravitational 
potential [36]. Interestingly, the same study found that a dis- 
crepancy did arise for high descent rates but not for high 
climb rates, descent being clearly more expensive than antici- 
pated. When flying steeply downwards at a descent angle of 
— 60°, the pigeons flew at a horizontal velocity of 3.6 m s -1 
and a vertical rate of descent exceeding 3 ms -1 . Hence, one 
possible explanation for this could be that flight costs 
during rapid descent are so minor and so rarely encountered 
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Figure 8. When the bird is ascending or descending, the quadrature phase 
arrangement is not always optimal. However, the phase 6 for which the power 
is minimized only deviates significantly from -90° during extreme rates of 
ascent or descent. The results here simulate a bird with a wingstroke frequency, 
f w = 7 Hz and equal amplitude oscillations on both forward and vertical axes 
of amplitude 8.5 mm. The locus of points at which the power is minimized at 
climb rates intermediate between those of the marked traces is also shown. 



error of quadrature phase model during ascent and descent 
40- 




n 1 1 r i i r 

-0.8 -0.6 -0.4 -0.2 0 0.2 0.4 
rate of ascent (ms _1 ) 



0.6 



o.: 



Figure 9. The quadrature phase assumption potentially overpredicts the body 
power for ascending or descending flight. The error is plotted here for climb 
rates within the range +0.9 m s~\ where the optimum phase appears to 
deviate linearly from 6= -90° in figure 8. The model is relatively accurate 
for sustainable rates of ascent but errors grow rapidly for high rates of descent. 



that there is little need for birds to acquire biomechanical 
flexibility that would allow significant deviations from quad- 
rature phase. However, this descent rate is more rapid than 
that which would be expected if the bird were simply gliding 
(typically no more than 2ms -1 for a pigeon). 

If more accurate estimates of biomechanical power are 
required then direct integration is an option, obtaining rela- 
tive velocities and displacements for the forward and 
vertical axes from which the individual energy terms can 
be derived, combining this information with altitude data 



(if available) to obtain body energy, and calculating power 
by numerically evaluating the mean rate of energy increase 
with time, being careful to ignore periods when the total 
energy of the bird is decreasing. This approach may offer 
improved precision, particularly when power is not all 
concentrated at the fundamental wingstroke frequency or 
when large departures from quadrature phase are expected. 



5. Discussion 

Despite the considerable complexities involved in flapping 
flight, by focusing on the consequences for the mechanical 
motions of the body, this oscillatory energy fluctuation 
model provides a useful initial step in theoretically underpin- 
ning the use of body-mounted accelerometers to estimate the 
relative costs of horizontal flapping flight in birds or bats. 
Acceleration-based proxies for the power detectable in the 
body of an animal during flight have been derived from 
first principles ((3.5), (3.7) and (3.9)) and in the appropriate 
units of ML 2 T~ 3 . Within these equations, the DBA formalism 
is encapsulated very naturally via r.m.s. acceleration ((3.14) 
and (3.15)), or aRMS. These expressions represent a substan- 
tial improvement over attempting to estimate the 
biomechanical energy expenditure or metabolic rate of 
flying animals from ODBA or VeDBA alone. The magnitude 
of ODBA varies with orientation, and different rescaling fac- 
tors must be applied to single axis projections of ODBA and 
VeDBA if they are to be used to estimate i/ rms or z rms (appen- 
dix B). We point out that aRMS, which is more closely related 
to VeDBA than ODBA, is just as straightforward to calculate. 
Furthermore, aRMS may be more universally applicable in 
future mathematical and empirical studies of animal loco- 
motion. DBA measures lack the units of power and so will 
always require a direct calibration against the rate of energy 
turnover. Their use in estimating the biomechanical costs of 
flapping flight has generally lacked a firm theoretical basis, 
pays no heed to sensor orientation relative to gravity, ascribes 
undue significance to accelerations in the horizontal plane 
and overlooks the significance of wingstroke frequency. As 
can be seen in (3.22), body power is rather sensitive at 
times to wingstroke frequency, so may in itself provide a 
valuable means of gauging flight effort independently of 
body power estimation. 

This analysis furnishes several novel parameters that 
may provide insights into the kinematics of flight. In par- 
ticular, it predicts that there may be two regimes of flight 
with regard to power production and wingstroke frequency 
and that there may be a transition from slow to moderate 
intensity flight, when equation (3.23) and figures 6 and 7 
suggest power savings may sometimes be possible. For 
steady horizontal flight involving sinusoidal vibrations on 
the forward and vertical axes, it has been shown that a 
quadrature phase arrangement is potentially advantageous. 
If flight strategies could exploit this, temporarily storing 
and retrieving energy during each wingbeat cycle, it might 
be possible for some species to avoid or curb what would 
otherwise be a third-order sensitivity in mean body power 
to wingstroke frequency. Due to the ubiquity of predators, 
the need to catch airborne prey and the general requirement 
for economy of locomotion, there may have been consider- 
able evolutionary pressure for birds to waste very little 
energy at wingstroke frequencies exceeding o) 0 . 



Nevertheless, body-mounted accelerometry is not privy to 
the subtleties of wing flexion, angle of attack and feathering, 
so this model leaves open the possibility that birds have 
considerable scope to adjust their flight style without 
necessarily compromising efficiency. 

When birds are flying in the linear regime then equation 
(3.14), which reflects the costs of combating gravity, shows 
that only the vertical, gravity-aligned component of the accel- 
eration should enter into the calculation of body power. 
Although forward accelerations become relevant to power esti- 
mation for the asymptotically cubic regime, equation (3.15) 
shows that their contribution to body power is subtractive. 
This somewhat counterintuitive need to subtract i/ rms from 
z rms in (3.7) stems from the fact that variations of gravitational 
energy do not always mask variations in kinetic energy. How- 
ever, quadrature oscillations in K y tend to erase fluctuations in 
Ky Z , not reinforce them. 

While the model directly considers vibrations only on the 
body of the bird, owing to mechanical coupling, the same 
kind of oscillations and trade-offs should also apply to the 
horizontal and kinetic components of the wings. However, 
because the centre of mass of the wings must travel a much 
greater distance during each wingstroke than the centre of 
mass of the body, for any given wingstroke frequency, 
fluctuations in wing kinetic energy grow quadratically with 
wingstroke excursion, but variations in wing gravitational 
potential grow only linearly. Therefore, the onset of the tran- 
sition between the linear and asymptotically cubic flight 
power regimes might be generally expected to occur at a 
lower frequency for the wings than for the body. The kin- 
ematics of the wings, which cannot be directly measured by 
the accelerometer, should generally dominate the overall bio- 
mechanical costs of flight. Indeed, kinetic energy fluctuations 
might easily become a more important consideration than 
compensating for the gravitational energy losses of each 
wingbeat. If the gravity terms are neglected, then (3.7) 
and (3.9) simplify somewhat and predict (P^) « 477 2 rab/ 3 
|B 2 — A 2 \. This could be especially true of energetic high- 
speed flight, even though there may then be significant 
wing retraction and supination to avoid undue aerodynamic 
drag [37,38]. Because the forces involved in wing retraction 
tend to cancel on the body due to bilateral symmetry, 
the costs are hidden from body-mounted accelerometers. 
Hence, the ratio of perceived body power to true total bio- 
mechanical power might be somewhat reduced during 
intense flight for some species, which may well require the 
tailoring of flight models to each species in the future, follow- 
ing empirical observations and extensions of the modelling. 

A number of other original summary statistics could 
prove useful to the interpretation of body-mounted accelero- 
metry data obtained from flying birds. The ratio A/B may 
help to characterize the mode of flight performance and per- 
haps also evaluate the skill and dexterity of individual birds. 
With the possible exception of hovering flight when birds 
may be able to recoup some of the kinetic energy stored in 
the air during the previous half-stroke, flying animals 
generally have no means of recovering energy lost to the 
environment. However, vibrations on the forward axis offer 
convenient energy storage which may also be exploited to 
reduce pitching of the body. Many birds use various reflex 
mechanisms during flight to subdue head vibrations and 
thereby avoid vision imp airment [33]. Body accelerations 
have a total amplitude Jy 2 + z 2 = a? VA 2 + B 2 which for 



small A/B can be approximated by o) 2 B(l + A 2 /IB 2 ). Hence, 
the additional loss of visual acuity due to a small A/B ratio 
would be relatively imperceptible, implying little need for 

A/B -> 0. 

The ratio A/B may also be particularly sensitive to effort in 
realistic situations, correlating with wingstroke frequency and 
increasing at higher forward velocities. The relative phase lag 6 
between the vertical and forward axes may act as a marker of 
ascending or descending flight, or reflect efforts to synchronize 
wingstroke frequency with other birds during V-formation 
flight. We also expect the ratios co/o) c , co/ co opt and 00/ co 0 to 
be informative regarding flight intensity and flight efficiency. 
Collectively, these measures may also crudely encode hints 
as to the altitude at which the bird is flying. Due to the com- 
plexities of wing kinematics and anatomical constraints, it is 
very likely that no simple unifying pattern will adequately 
summarize all species, but departures from normality are 
often the most interesting aspects of biological research and 
so additional parameters can prove very useful in highlighting 
departures from non-conformity. Therefore, these flight vari- 
ables may be particularly valuable in helping to unpick the 
challenges involved in flying efficiently. The static acceleration 
also offers a potentially illuminating variable for flying ani- 
mals which has been largely ignored to date. Birds can 
sustain prolonged banking when circling or jostling for pos- 
ition within a cluster flock [39], and any drift in the mean 
direction of the momentum vector induces a non-gravitational 
contribution to the static acceleration. Thus, when the static 
acceleration deviates appreciably from gravity, it would 
suggest that the bird is not undergoing steady horizontal 
flight. However, the converse is not true, because one also 
expects the static acceleration to tally with gravity during 
steady non-banking ascending or descending flight. Therefore, 
the distribution and time variability of the static acceleration 
can be informative. 

Dimensional considerations may allow the results 
obtained here to be extrapolated to some degree, parti- 
cularly regarding the estimation of biomechanical power 
from accelerometry for aquatic animals. Due to the buoy- 
ancy afforded by water, the estimation of inertial costs for 
aquatic species during swimming is not encumbered by 
gravitational considerations [2]. This invulnerability to grav- 
ity suggests that the cost of swimming should correlate with 
the product of body mass, the period of the swimming 
stroke and some function of the square of the decomposed 
r.m.s. accelerations, dependent upon the anatomy of the 
species under consideration. Locomotion costs in terrestrial 
animals are likely to be more complex: weight support 
can either be provided continuously or episodically by 
the ground. 

Naturally, there are limitations to what a body-mounted 
accelerometer alone can glean about flight. During free- 
ranging flights, there could be circumstances where basic 
inferences may be misleading, particularly if the rate of 
ascent or descent is unknown. In addition, special care may 
be needed when attempting to disentangle the static and 
dynamic accelerations for birds using intermittent modes of 
flight such as flap-gliding or flap-bounding. Accelerometers 
cannot infer absolute velocities in any direction, and many 
birds exploit the assistance of thermals, following winds 
and airflow over uneven terrain, all of which are capable 
of drastically altering the power requirements of flight. 
Nevertheless, the present model offers a practical and non- 



invasive method of extracting from accelerometry a variety 
of parameters that could be informative concerning flight 
style and performance, while also providing an explicit 
procedure for determining biomechanical body power in 
free-flying birds which may be generally proportional to 
overall flight power. In the complex processes that transform 
the biochemical energy of birds into atmospheric vortices, 
turbulence and heat, aerodynamic costs lie downstream of 
the biomechanical costs. While it might eventually be poss- 
ible to incorporate them within an extended model, the 
formidable challenges of contending with complicated 
wake patterns, vortex interactions and chaotic flow patterns 
continue to plague theoretical models, and quasi-static 
approximations to the Navier- Stokes equations commonly 
used in aerodynamic analysis inherently limit their accuracy 
and usefulness [40]. A more realistic near-term goal would 
be to broaden the present modelling to include wing kin- 
ematics and morphology. This will inevitably necessitate 
the input of anatomical information allowing the body 
power relationship to accommodate allometric differences 
between species and also address wing-propelled locomotion 
in the media of differing densities. Due to the intrinsic com- 
plexities, we anticipate that experimental data collected from 
a variety of species will be required. Birds are graceful aero- 
nauts, skilfully adjusting their posture and technique in 
ways we have only started to perceive [41]. However, a 
battery of physiological, biomechanical and aerodynamic 
techniques can augment and refine one another when quan- 
tifying flight costs. These complementary approaches to 
studying avian energetics hold much promise in arriving at 
a more unified understanding of the compromises involved 
when animals fly — whether they are foraging for food, 
migrating, chasing airborne quarry or evading predators. 
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Appendix A. Reorientation during 
post-processing 

When flying in a straight line, the static acceleration 
a s = (#1,^2,^3), where ak represents the time-averaged 
mean of the kth accelerometer channel, should equal the 
gravitational acceleration, g. When this does not hold, it is 
possible that the bird is turning, undulating, bounding or 
flying through zones containing vertical air currents. Centri- 
petal acceleration combines with gravity during turns, the 
direction of the resultant vector determining the degree of 
banking necessary in order that the bird experiences the 
net acceleration dorsally so that the forces on the wings 
are symmetrically balanced. At such times, the increase in 
the static acceleration relative to gravity is a useful guide- 
line as to the departure from linear motion. For uniform 
rectilinear motion, the mode of the distribution of the 
static acceleration magnitude can be identified with gravity. 
To some degree, this allows for the self-calibration of accel- 
erometry data. 

When ||a s || ~ ||g||, a condition that is easily checked, it is 
straightforward to calculate the vertically aligned com- 
ponent of the acceleration a vert using the scalar projection 



of a onto a q : 



fli fli + a 2 a 2 + a 3 a 3 a\U\ + a 2 a 2 + 03^3 , A - x 

flvert = J =— ~ • (Al) 



^Ja\ + a 2 + a; 



This vertical acceleration is the combination of both static 
and dynamic components, so the dynamic vertical accelera- 
tion a z can be obtained by subtracting the time-averaged 
value of a vert using 



Q-z — ^vert ^vert- 



(A 2) 



Projections of a orthogonal to the vertical axis will then lie 
in the horizontal plane: 



and 



a\a 2 + a 2 a 3 + a$ai 



a\a 3 + a 2 a\ + a$a 2 



(A3) 



(A 4) 



As before, the dynamic acceleration is obtained by 
subtracting the static acceleration (#hi = cihon — #hori and 
#h2 = #hor 2 —^hor 2 )' A means of determining # x and a y from 
«hi an d «h2 is then required. By symmetry, one expects the 
dynamic acceleration of the forward (y) axis to exceed that 
on the lateral (x) axis. One way to proceed would be to 
first determine the fundamental wingstroke frequency 
using a z then bandpass filter a hl and a h2 using fast Fourier 
transforms so that only frequencies near the fundamental 
wingstroke frequency are retained. Following this, one 
could compute the angles cf) = arctan(ah2/^hi)/ compile a 
circular histogram of the results and use the angle corre- 
sponding to the peak in the histogram <P, to reorient a hl 
and # h 2 as follows: 



and 



a x = «hi si n @ — a h2 cos <P 



(A 5) 



(A 6) 



There is a potential ambiguity in this result concerning 
the polarities of a x and a y . It arises due to the fact that 
one expects two peaks in the angular histogram separated 
by 77. Therefore, one should also evaluate a x and a y using 
0 <— 0 + 77. In practice, there may be no need to calculate 
a x because it is unlikely to represent interesting information, 
but the polarity of a y is potentially important if one is eager 
to know, for example, how the phase shift between the for- 
ward and vertical axes varies with time. This model expects 
that the displacement on the y-axis will always lag behind 
that of the z-axis, and hence a z should always lead a y , resol- 
ving the ambiguity. One can again apply bandpass filtering 
to a y and a z around the detected wingstroke frequency in 
order to test which value of <P is appropriate. The phase 
shift can be accurately measured in the recovered time 
domain after bandpass filtering using linear interpolation 
between samples in the vicinity of the positive /negative 
going zero crossings. If the orientation of the accelerometer 
is fixed with respect to the body of the bird, then this pro- 
cess need only be performed once, and the value of <P can 
then be reused without recalculation. 



Appendix B. Relationship between overall 
dynamic body acceleration, vectorial dynamic 
body acceleration and r.m.s. acceleration 

Consider the instantaneous dynamic acceleration a d = (a x , Oy, a z ) 
derived from a triaxial accelerometer. The generalized mean or 
iAnorm of the components of a& is defined as 

\v\Vv 



+ to/ r + \<h 



(Bl) 



Therefore, the magnitude of according to VeDBA is 
a x + a l + a h whereas its magnitude according to ODBA is 



\&x\ + \ a y\ + \ a z\- Converting to spherical coordinates (r, 6, (/>), 
one has 



rsin 0cos <fi, 
■ r sin 0 sin 4> 
- r cos 0. 



and a z 

The magnitude of a& using VeDBA is thus 

" 2 = (r 2 sin 2 0cos 2 (/> + r 2 sin 2 0sin 2 4> 



(B2) 
(B3) 
(B4) 



at + at + at 



+ r 2 cos 2 
= r. 



W2 



(B5) 



This confirms the standard expectation of Euclidean 
trigonometry and the Pythagorean theorem. However, the 
magnitude of a& in the case of ODBA generally disagrees 
with this because 



and 



/ ODBA V 
\ VeDBA J 



/ODBA V 



a x + \a v 



al + a 2 + al 



V VeDBA/ 



1 + 2 



gxgyl + Igygzl + |gA| 

al + al + al 



(B6) 



(B7) 



Due to the second term, ODBA will exceed VeDBA unless 
at least two of the three acceleration components are zero. 
The two will agree only when there is alignment of the accel- 
eration with one of the three measurement axes. Hence, the 
response of ODBA varies according to orientation. 

To find the maximum error in ODBA, let a y = a x + a and 
a z = a x + p where, without loss of generality, it can be 
assumed that neither a x/ a y nor a z are negative. Now let £ = 
(ODBA 2 / VeDBA 2 - l)/2 so that the stationary points of £ 
will be identically located to those of ODBA/ VeDBA: 

a x a y + a y a z + a x a z 3a x + 2aa x + 2pa y + ap 



3a x + 2x(a + 0) + a 2 + f3 2 



(B8) 



Equating to zero the partial derivatives of £ with respect 
to a and /3 gives 



(a 2 + a 2 + a 2 z )(2a x + /3) = 2(a x + a)(a x a y + a y a z + a x a z 



and 

0 2 + a 2 + a\){2a x + a) = 2(a x + P)(a x a y + a y a z + a x a z 

These expressions reduce to 

(a x + a)(2a x + a) = (a x + p)(2a x + j8) 

and 

(a + 0)(a-p) = -3a x (a-0). 



(B9) 

(BIO) 

(Bll) 
(B12) 



The first condition is satisfied only if a = f3, which also 
ensures the second is satisfied. When a ^ fi, the second 
condition holds if a + /3 = — 3a x which cannot be true given 
that neither a y nor a z axe negative, which implies that 
a + p > -2a x . This leaves only a = /3 and hence \a y \ = \a z \ 
after restoring moduli. For \a x \ —> 0, one then finds that 
(ODBA/VeDBA) 2 can be no less than 2. This can be seen 
by letting y = \a y \ = \a z \ and \a x \ = e with £<y: 

\ + 2( tt l€ t \ ^2 + 2e/y>2. (B13) 



/ ODB A \ 
\ VeDBA J 



2y 2 + e 2 



However, as \a x \ —> \Oy\ = \a z \ one finds that (ODB A/ 
VeDBA) 2 can be no greater than 3. This can be seen by 
letting \a x \ = y ± e where, once again, e<y and y = \dy\ = \a z \: 



/ ODBA V 
\ VeDBA J 



1+2 



1 + 



3T 2 ± ley 
3T 2 ± 2ey+e 2 



1 + 673T 2 



<3. 



(B14) 



Hence, ODBA/VeDBA is maximized when the accelera- 
tion components on each axis have an identical magnitude. 
Because the minimum is already known to occur when two 
of the components are zero, one may conclude that ODBA 
is confined to the range 

VeDBA < ODBA < VeDBA. (B15) 

It is interesting to ask whether VeDBA can be recovered 
from historical records of ODBA. This is possible only in 
an approximate statistical sense, and best results would 
be obtained when data have been collected from animals 
whose orientation in space varies considerably or in situ- 
ations where the orientation of the accelerometer itself is 
free to drift. In order to determine the rescaling factor, it is 
first necessary to calculate the mean exaggeration of ODBA 
relative to the true acceleration magnitude. Due to symmetry, 
it is sufficient to consider the solid angle ft = tt/2 corre- 
sponding to the octant 0 < 0 < it 12, 0 < <f> < tt/2. Although 
the mean value of VeDBA is simply the vector length r, the 
mean value of ODBA is 



(ODBA) 



2 -\ 



ODBA dil : 

■/2 ,tt/2 



[\a x \ + \a y \ + \a z \)df2, (B16) 



2 PTT/Z PIT/ 

(rsin 0cos 4> + rsin 0sm4> 

'""Jo Jo 

') sin#d(/>d/9, 
[sin 2 #sin 4> — sin 2 0cos cf> 



2r 

77 



rcos ( 

tt/2 



tt/2 



1 



2r 

77 

2r 

77 



2r /77 77 

: ^(2 + 8 + 



cos 20 + — sin 20 
4 



d6, 



sin 20 77 cos 20 



»W 2 



3r 
~2 ' 



(B17) 



(B18) 



(B 19) 



(B20) 



Therefore, on average, ODBA exaggerates the true accel- 
eration magnitude by 50% and so, in some circumstances, 
ODBA data records can be translated into estimates of 
VeDBA simply using 

VeDBA « \ ODBA. (B21) 



A highly significant linear relationship between these two 
measures has already been experimentally observed, with the 
best fit corresponding to VeDBA w 0.014 + 0.6418 ODBA for 
units of g [27]. The slope of this empirical relationship agrees 
with the theoretical value to within 4%. In the same work, the 
envelope of figure 2 exhibits a wedge distribution whose 
upper and lower slopes are approximately unity and l/\/3, 
corresponding to the anticipated range in error of ODBA 
due to changes in orientation. 

When evaluating r.m.s. accelerations, there is no explicit 
requirement to calculate vector lengths, but we stress that 
the Euclidean formulation, as adopted by VeDBA, is implicit 
in the present modelling. While VeDBA and r.m.s. accelera- 
tion are in agreement here, and one would generally expect 
to find an excellent correlation between them, their magni- 
tudes are anticipated to differ whenever there is any spread 
in the distribution of the dynamic acceleration data. For com- 
plex or aperiodic acceleration profiles, we recommend the 
r.m.s. method over VeDBA in all cases, because the am- 
biguity implies that no single rescaling factor will suffice. 
However, for sinusoidal motion if/= f3 sin(f ) along a straight 
line, a simple rescaling is possible. The time-averaging used 
by VeDBA follows that of ODBA, namely the I^-norm or 
arithmetic mean. However, that of the r.m.s. prescription 
follows the L 2 -norm or quadratic mean which, as found 
previously in (3.11), yields (iA rms ) = P/\/2. If follows 
from the properties of the generalized means that (t/fyeDBA) 
will never be smaller than (*A rms ) and its value is 



VeDBA 



sin(f)df = [ 

0 77 



-cos(0]5 



2/3 

77 



(B22) 



For this simplified situation, one finds 2v / 2(iA rms ) ~ 
^^/VeDBA) and hence 

aRMS « VeDBA « -^L ODBA. (B 23) 

In applications where approximate alignment of the accel- 
erometer's z-axis with gravity is attempted, small changes in 
orientation occurring at frequencies that are not rejected by 
the sampling window will cause ODBA to increase even in 
the absence of vibration. Because the magnitudes of a x and 
a y are then much smaller than gravity one finds that 



/ ODBA y 
V VeDBA; 



1 | 2 lgx|_+|gy| 



(B24) 



Hence, the lack of rotational invariance in ODBA (as 
depicted in figure 10) then causes it to respond linearly to tilt- 
ing of the accelerometer, tending to accentuate the sensitivity 
of ODBA to activity. Dependent on the moment of inertia 
there can be a cost associated with periodic adjustments in 
attitude, as might occur in animals exercising on a treadmill, 
and ODBA may be sensitive to it. For arbitrary orientations 
with respect to gravity, ODBA can both exaggerate and 
underestimate changes in the acceleration magnitude 
caused by rotation. Almost any acceleration measure is vul- 
nerable to errors when the rotation of a transducer is not 
limited to low frequencies because it is then impossible to 
accurately separate the static and dynamic accelerations with- 
out information from a gyroscope. 

When significant rotation exists, in the absence of a 
gyroscope, some measure of the variability of the instan- 
taneous magnitude of the total acceleration vector a could 



to 
b 



CD 

CD 
-1^ 




Figure 10. Variation of the error in ODBA with orientation, rising from a minimum of zero (six centres of the square contours) to a maximum of 73.2% (eight 
centres of the circular contours) with a mean error of 50% (contours at 10% intervals). A single contiguous contour is obtained when ODBA/VeDBA = y/i. (Online 
version in colour.) 



offer an alternative proxy for biomechanical effort. For 
in stance / one m ight calculate the standard deviation of 
^jAl +A*+Al for all the raw outputs of the accelerometer 
falling within some time interval. This could prove 
informative whenever a significant component of animal 
activity involves body rotation. Because accelerometers 
remain sensitive to gravity even within buoyant media, 
this may be especially useful in the context of aquatic 
animals as they are not constrained by gravity when swim- 
ming underwater. In circumstances where the dynamic 
acceleration is extremely small, activity can also be 
estimated by quantifying the rate of body reorientation 
d4>/dt. This can be approximated as Arf>/ At using pairs 



of static acceleration vectors A s \ = ( A x \ , A y \ , A z \ ) and 
A s2 = (A x2 ,A y2 ,A z2 ) separated by a fixed time interval At 
appropriate to the rotation rates of interest: 



to 
b 



CD 

CD 
■4^ 



A<f) 



IjAsi x A s2 || 
[ \\A s i\\-\\A s2 \ 



(B25) 



and 



Ac)) = arcsin a 



(A yl A z2 - A y2 A zl f + (A xl A z2 - A x2 A zl f 
+(AxiA y2 - A x2 A y 



(B26) 
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